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The model of Fermi particles with random two-body interaction is investigated. This model allows to study the origin and 
accuracy of statistical laws in few-body systems, the role of interaction and chaos in thermalization, Fermi-Dirac distribution 
for quasi-particles with spreading widths, matrix elements of external field and enhancement of weak perturbation in chaotic 

compound states. 
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As is known, excited states in many-body systems such 
as compound nuclei, rare-earth atoms, molecules, atomic 
clusters, quantum dots in solids, etc. are very compli- 
cated ("chaotic") and can be described via statistical 
methods. There are two major approaches. The first 
is based on random matrix theory (see e.g. JjJ). This 
approach is very general in nature and therefore it does 
not describe many important properties of realistic many 
body systems, its prediction being limited to level statis- 
tics, localization properties of eigenstates and the like. 
The other approach is based on direct numerical investi- 
gations of the given particular system. For example, the 
results of direct diagonalization of the Hamiltonian ma- 
trix for the rare-earth Ce-atom have been compared with 
statistical theory of compound states [^]. A similar study 
has been performed for the s — d shell model of complex 
nucleus |3| where the problem of thermalization has also 
been considered. Quite obviously with this second ap- 
proach, it is very hard (e.g. due to the lack of statistics) 
to draw general conclusions concerning the accuracy of 
statistical laws in systems with finite number of particles 
and the conditions for their applicability. 

Here we suggest an " intermediate" approach based on 
a simple mathematical model with random interaction, 
which however takes into account the most important 
features of many-body systems such as single-particle 
orbitals, two-body interaction and the Pauli principle. 
In our model of Random Two-Body Interaction (RTBI) 
we assume that the system consists of n Fermi-particles 
which can occupy m orbitals with no double occupancy 
(extension to Bose-systems is straightforward). As an 
example, the energies of orbitals are taken in the form 
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where the second term is introduced in order to avoid 
the degeneracy of many-particle states. Matrix elements 
of two-body interaction are chosen as random variables 
distributed according to the Gaussian law with the zero 



mean and the variance Vq . Therefore, our model has 3 
independent parameters: n, m and V$/d . 

Structure of the Hamiltonian matrix. The basis 
of our Hamiltonian matrix is chosen to correspond to the 
products of single-particle orbitals ordered by increas- 
ing many-body unperturbed energies. Therefore, with- 
out interaction the matrix is diagonal, with increasing 
elements. The size N of RTBI-matrix is given by N = 
C™ = ml/(n\(m—n)\) and is exponentially large for large 
n and m . Due to the two-body character of interaction 
the matrix elements are zero between those basis states 
which differ by more than two orbitals. Therefore, the 
Hamiltonian matrix is essentially sparse with the sparsity 
(ratio of non-zero elements to the total number N 2 ) given 
by s = (1 + n(m — n) + n(n — l)(m — n)(m — n— l)/4)/7V 
At large n, m the sparsity is exponentially small, 



s ~ exp(— n) . In fact, the number N% 
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of independent variables (number of pair interactions) in 
the model is even smaller than the number of non-zero 
Hamiltonian matrix elements. 

There are three types of interactions between many- 
body states: diagonal interaction which contains ki nt = 
n(n — l)/2 two-body terms, interaction between states 
which differ by one orbital with ki n t = n — 1 terms and 
interaction between states which differ by two orbitals 
with ki n t = 1 ■ Finally, distant many-body states which 
differ by more than two orbitals have zero matrix ele- 
ments. Thus, the Hamiltonian matrix has large and 
increasing diagonal, plus sparse band-like structure with 
a decrease of off-diagonal elements as a function of the 
distance \i — j\ from the diagonal. 

In what follows we consider a particular case which can 
be used to describe Ce-atom : m = 11, n = 4 (therefore, 
N = 330), d = 1,Vq = 0.12 (the two latter parameters 
are given in eV). Direct investigation of this atom (see 
details in ^|) has shown that it can be treated as a very 
chaotic system. 

Energy spectrum and eigenstates. As was found 
in 0] (see also ||), for two-body random interaction the 
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density of states (DOS) should be of Gaussian type, pro- 
vided m 3> n 3> 1. Our numerical data show that in 
the RTBI model in spite of small number of particles 
the DOS is quite close to a Gaussian with (E) = 25.1 
and a = 5.7. Another commonly discussed characteristic 
of the spectrum is the level spacing distribution; in our 
model it is well described by the famous Wigner-Dyson 
law. This fact is also in agreement with old studies ||] 
of some two-body interaction models (see also references 
in H). However, our interest goes beyond the DOS and 
spectrum statistics. 

A quantity of relevant interest is the so-called spectral 
local density of states or strength function. It is defined 
by the weight w(E, Ej) = |C?| of a particular basis com- 
ponent j = 1,2, ..,7V in the eigenstates with an energy 
close to E , according to the relation 



p(j,E) = 



w(E,E J 
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Here D is the local mean level spacing. The width of 
this function in energy E is proportional to an effec- 
tive number of components in the expansion of an un- 
perturbed basis state in terms of exact eigenstates. In 
nuclear physics, this function is usually assumed to have 
the Breit-Wigner form, p(j, E) ~ [(E - E,j) 2 + T 2 /A] _1 
(the relevance of this function to ergodicity and chaos is 
studied in According to our numerical results, in the 
center of the spectrum (j w N/2) the spreading width is 
equal to T c ss 1.0. One should note that close to the 
edges j = 1, N , the symmetric form of the distribution 
p(j, E) is strongly distorted. However, the width of the 
distribution itself changes slightly. What is more impor- 
tant, the tails of the function p(j,E) decay much faster 
than in the Breit-Wigner law. This is the consequence of 
the band-like structure of the Hamiltonian matrix. Both 
the value of F c and the fast decay of the tails, are in 
agreement with Ce-atom calculations Q . 

The localization length (number of principal compo- 
nents of the eigenstates) can be defined through the re- 
lation I ~ exp ((S 1 )) where S is the statistical entropy of 
individual eigenstates (see, e.g. ||]). One should note 
that in spite of a completely random character of the 
interaction, the lowest states turn out to be quite sim- 
ple containing I « 1 basis components. This can be 
explained by the low density of states near the ground 
state. The localization length I is maximal in the cen- 
ter E w (E) of the energy band and, according to our 
numerical data, can be well described by the Gaussian 
function: I = Acxp(-(E - {E)) 2 /(2a 2 )) with A w 135 
and <7 w 5.45 .Thus, in the center of the spectrum the 
number of principle components is about 100 which is 
again very close to the results of direct calculations for 
Ce-atom @. 

Statistical treatment of finite Fermi-systems. 

As is known, quantum statistical laws are derived for sys- 
tems with infinite number of degrees of freedom, or for 
systems in a thermostat. From this point of view, it is of 



importance to study how statistical laws appear in sys- 
tems with finite number of particles. Below, we show that 
in the RTBI model one can introduce a reasonably ac- 
curate statistical description based on such macroscopic 
characteristics as the temperature T, the chemical po- 
tential p and average occupations numbers rii for the or- 
bitals. In the mean field approximation these quantities 
can be obtained from the following set of equations: 
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n(e a ) = (1 + exp(e Q - p)/Ty 



< ^ Va P n < 3 > 
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In our case, the average interaction < V a p > is zero, 
therefore, we can omit interaction term in (|^) and (|^). 
On the other hand, there exist other important effects of 
the interaction which appear beyond the mean field ap- 
proximation. Namely, the interaction leads to the spread- 
ing width for the basis states and for the single-particle 
orbitals (r Q ) . It also results in the shift of average ener- 
gies, e a — e a + 6e a . According to our numerical data, 
the magnitudes Se a are smaller than r a and vanish in 
the mean (Se a < for low orbitals, 8e a > for high or- 
bitals and Se a « near the center). For this reason, we 
will take into account the effect of the spreading width 
T Q only. 

Instead of (^), by averaging the Fermi-Dirac distribu- 
tion n(e) over the interval r a : 



e a +r a /2 
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we now introduce the occupation numbers (g) which 
take into account the finite spreading width of "quasi- 
particles" . In the limit r a = the expression (|J) with 
n a = n(e a ) is recovered. The numerical solution of Eqs. 
([|||) is presented in Fig.l. In order to reveal the influ- 
ence of the spreading width r Q , two curves are given for 
comparison with r a — and T Q = 3.0 . The latter value 
of r a was deliberately taken larger than the spreading 
width of basic components T c « 1.0 in order to elucidate 
the effect of r Q . The temperature T and chemical po- 
tential p are calculated from the above equations (j^Jg) 
at fixed excitation energy SE = E — E m i n . We found 
that for small 5E the temperature is quite different for 
different values of r Q . Surprisingly, in spite of this, the 
two "theoretical" curves practically coincide. This means 
that the temperature mimics the statistical effect of the 
interaction, the phenomenon far from being trivial. 
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The average occupation numbers directly computed 
from exact eigenstates are shown in Fig.l with circles. 
One can see that even for 4 particles the actual distri- 
bution of occupation numbers can be approximately de- 
scribed by the statistical Fermi-Dirac distribution. How- 
ever, there is a clear deviation which indicates that for 
the chosen parameters the thermalization is not com- 
plete. Numerical data show that the deviation disap- 
pears with increasing the excitation energy SE , while it 
increases when the perturbation Vb decreases. The latter 
effect has been also observed in the s — d nuclear shell 
model || where the relevance of chaos to the thermal- 
ization was studied in a different approach. We have to 
stress that equilibrium distribution, or "thermalization", 
in this few-particle system is due to the interaction which 
"chaotically" mixes neighboring basis states with differ- 
ent occupation numbers. To demonstrate this, we plot in 
Fig. 1 the occupations numbers for the same system with 
no interaction (Vo = 0) . Their distribution has nothing 
to do with the Fermi-Dirac distribution: it is singular 
and for some values of SE even not monotonous. 

Correlations between occupation numbers. In 
a few body-system one could expect quite strong cor- 
relations between occupation numbers of different or- 
bitals. We have found, instead, that typically the cor- 
relations are weak; even for close orbitals the ratio 
(h a ) (ftp) I (n a np) ks 0.82. Only when occupation num- 
bers are small, the correlations can be very strong. 

Matrix elements of an external perturbation. 
The main problem in the compound state theory is the 
calculation of effects of an external perturbation. Since 
matrix elements of any single-particle operator can be ex- 
pressed in terms of elementary density matrix operators 
Pa/3 — O'aO'p which transfer the particle from the orbital 
to the orbital a , the main interest is in statistical prop- 
erties of matrix elements p a p . Below, we use recent ap- 
proach developed in |^] (see also the study of Ce-atom in 
where the following expression for the mean square 



matrix elements (MSME) ^ ni) = 
compound states (ni|and (n 2 |has been derived: 



(n 2 ni 
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where w Q/3 = e a - tp and Q aj3 = {n 2 \n a (l ~ h/3)\n 2 ). 
Here the sum runs over many-body basis components 
with w ni , w n2 being the weights of these components in 
the states n 1 and n 2 (see Eq. (||)). If we know the shape 
of the strength function w, this sum can be replaced by an 
integral and explicitly evaluated. In particular, when w 
iis Breit-Wigner type, the final expression for the MSME 
has also the Breit-Wigner form (see details in J2|). 

In order to check the accuracy of the above expres- 
sion (R) in the RTBI model, we have directly calculated 



f° r a ~ 4; /3 — 5 (transition from the ground 
state to the nearest one) and different values of m and 
n 2 . To reduce fluctuations, an averaging over a num- 



ber of realizations of the Hamiltonian matrix has been 
done. The direct comparison with the analytical pre- 
diction (0) shows a quite good agreement. In particu- 
lar, the positions of maximum as well as the widths of 
^ofl 2 ™ 1 ^ 01 dependence on n\ and n 2 are well described 
by (f7|). However, quite unexpected deviations have been 
discovered which were found to be generic in the model. 
A typical example is given in Fig. 2 where a clear differ- 
ence is seen between statistical approximation (0) and 
numerical data. The absolute difference is maximal at 
the center of the ^-dependence; however, in the tails the 
relative difference is even larger. The thoroughly study 
of this phenomenon (see details in [|l0||) has revealed a 
very intriguing fact: the origin of this effect is in the 
underlying correlations induced by the two-body nature 
of interaction. Similar correlations were observed in the 
model of random separable interaction [jll],[l2) ■ Full ana- 
lytical treatment of the correlations for the RTBI model 
is given in Jl0| . In particular, for the tails the contri- 
bution of the correlation term £ corr to the total MSME, 
t;totai = £,corr + £stat , has been estimated as 



I? 



£,stat 



(n — 2)(m — n — l)(m — n + 2) 
n(m — n)(m — n + 3) 



(8) 



where £, s tat is very close to that given by the expression 
(^). In the maximum of £, the estimate for m - n » 1 
reads as £,totai / istat ~ 2-2m/(n(m-n)). The remarkable 
result is that these correlations do not decrease with an 
increase of number of particles n and number of orbitals 
to. Numerical data for larger values of n — 7 and to = 14 
(with N — 3432) have confirmed this prediction. 

Enhancement of a weak perturbation in chaotic 
many-body systems. The RTBI model allows for the 
study of a very important effect, namely, the enhance- 
ment of a weak perturbation W, which takes place in 
systems with chaotic compound states. This effect is pro- 
portional to the mixing coefficient rj —< ni\W\n 2 > / Ai 2 
where A12 = E± — E 2 is the spacing between the neigh- 
boring energy levels of compound states. Since the spac- 
ings for strongly excited states of many-body systems 
are exponentially small, A12 ~ exp(— n) , one could ex- 
pect strong enhancement of the perturbation in compar- 
ison with the single-particle mixing defined by rj s =< 
a I W \(3 > /(e a — £/?)■ The possibility of the enhancement 
in compound nuclei was pointed out in [ p^|p^ | and consid- 
ered in details in recent review However, such non- 
trivial effects like repulsion between energy levels and the 
above discussed correlations, may have strong influence 
on the enhancement. Our preliminary results show that 
even for 4 particles the enhancement in the RTBI model 
does exist. 

Conclusive remarks. In this Letter we have anal- 
ysed the RTBI model which, unlike conventional ran- 
dom matrix models, allows for the study of many im- 
portant problems related to the two-body character of 
the interaction. As an example, the role of interaction 
for the appearance of the Fermi-Dirac distribution has 
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been investigated for the parameters of Ce atom. It was 
found that two-body interaction gives rise to thermaliza- 
tion and that the statistical effect of interaction can be 
imitated by an increase of temperature. The study of cor- 
relations between occupation numbers has revealed suffi- 
ciently weak correlations even for small number (n = 4) 
of particles. This justifies the approximation of indepen- 
dent particles which is typically used in the description 
of complex compound states. To describe the effect of 
spreading widths of orbitals, the generalization of the 
Fermi-Dirac distribution was suggested. 

The numerical and analytical treatment has shown 
that the statistical theory reproduces quite well the 
global structure of matrix elements of an external per- 
turbation between compound states. On the other hand, 
underlying correlations have been discovered which are 
induced by the two-body character of interaction, even 
if the latter is completely random. This phenomenon 
results in serious deviations from the statistical predic- 
tions for matrix elements between compound states. At 
the moment, all consequences of these correlations are 
still not understood, however, preliminary data jl5| ] show 
that they can lead to the so-called gross-structure (sharp 
peaks) in a cross section. 

In conclusion, our results show that the discussed 
above TBRI model can be a very useful tool in the study 
of many important problems of statistical physics of com- 
plex quantum systems. 
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FIG.l. Distribution of occupation numbers for the exci- 
tation energy 6E = 1.63 averaged over an ensemble of 5 
Hamiltonian matrices. The full and dotted lines are nu- 
merical solutions of Eqs. (|]]|) with T a = and T a = 3 
respectively. The circles and squares give results of the 
direct computation of average occupation numbers with 
and without interaction respectively. 

FIG. 2. Mean square matrix element £ calculated for 
rii = 55, a = 4, f3 = 5 as a function of n<i. Averaging 
over 100 Hamiltonian matrices with different realization 
of random two-body interaction has been made. The dots 
correspond to the direct numerical computation, the solid 
line represents the statistical approximation (R) . 
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